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Abstract 

We use the quantum Brownian model to derive the uncertainty relation for a quantum 
open system in an arbitrarily-squeezed initial state interacting with an environment at 
finite temperature. We examine the relative importance of the quantum and thermal 
fluctuations in the evolution of the system towards equilibrium with the aim of clarifying 
the meaning of quantum, classical and thermal. We show that upon contact with the 
bath the system evolves from a quantum-dominated state to a thermal-dominated state 
in a time which is the same as the decoherence time calculated before in the context of 
quantum to classical transitions. We also use these results to deduce the conditions when 
the two basic postulates of quantum statistical mechanics become valid. 
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It is a well-known fact in quantum mechanics that a lower bound exists in the product 
of the variances of pairs of noncommutative observables. Taking the coordinate x and 
momentum p as examples, the Heisenberg uncertainty principle states that with (Ax) 2 = < 
x 2 > — < x > 2 , and (Ap) 2 = < p 2 > — < p > 2 , the uncertainty function 

h 2 

Uq M = (Ax) 2 (Ap) 2 > — (T = 0, quantum mechanics) (1) 

The existence of quantum fluctuations is a verified basic physical phenomenon. The origin 
of the uncertainty relation can be attributed as a mathematical property of Fourier analysis 
[1] which describes quantum mechanics as a wave theory. Recent years have seen effort in 
establishing a stronger relation based on information-theoretical considerations [2] [3] . 

In realistic conditions quantum systems are often prepared and studied at finite tem- 
peratures where thermal fluctuations permeate. At high temperatures the equipartition 
theorem of classical statistical mechanics imparts for each degree of freedom an uncer- 
tainty of kT/2. Thus the uncertainty function for a one-dimensional harmonic oscillator 
approaches the limit 

kT 

U? 3 ~ (tz~ ) 2 (high T, classical statistical mechanics) (2) 

where hQo is the energy of a normal mode with natural frequency fio- This result, ob- 
tained by assuming that the system obeys the Maxwell-Boltzmann distribution, is usually 
regarded as the classical limit. For a system of bosons in equilibrium at temperature T, 
the application of canonical ensemble gives the result in quantum statistical mechanics as 

U^ E = — [coth(— — )] 2 (all T, quantum statistical mechanics) (3) 

which interpolates between the two results (1) and (2) at T = and T » hflo/k. This 
result applies to a system already in equilibrium at temperature T. 

Our purpose here is to study the corresponding non-equilibrium problem. At time to 
we put the system in contact with a heat bath at temperature T and follow its time evolu- 
tion. We want to see how the uncertainty function L/r(£) changes from the initial quantum 
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fluctuation-dominated condition to a later thermal fluctuation-dominated condition. By 
comparing this result with the decoherence studies recently carried out [4,5,6], where two 
characteristic times — the decoherence time tdec and the relaxation time t re i — are defined, 
one can apply the physics of these two processes involved to examine the following three 
interrelated issues : 

1) the realization of the basic tenets of quantum statistical mechanics from quantum dy- 
namics; 2) the relation between quantum and thermal fluctuations and 3) their role in the 
quantum to classical transition. 

Quantum statistical mechanics of a macroscopic system is derived from the quantum 
dynamics of its microscopic constituents under two basic postulates [7]: i) random phase, 
and ii) equal a priori probability. The first condition enables one to assign probability 
distributions to a system occupying certain quantum states. It requires sufficient dimin- 
ishing of interference terms in the wave function or that the density matrix of the system 
be approximately diagonal. The second condition when applied to a closed system forms 
the basis of the micro canonical ensemble. When applied to an open system in contact with 
a heat bath described by the canonical ensemble, it ensures that the system equilibrates 
with its environment. We want to examine the processes by which these two conditions 
are attained from a more primitive level, starting with the microdynamics of a system of 
quantum particles. Specifically, we want to see if there is a characteristic time when the 
phase information is lost (Postulate i) and another time when the system attains equi- 
librium with its surrounding so that all accesible states in the closed system are equally 
probable (Postulate ii). 

On the second issue, the demarkation of classical, quantum and thermal regimes are 
not always clearly noted, their usual definitions or usage oftentimes are imprecise. With 
the aid of the results obtained here we will point out some existing confusions and make 
their meanings more precise. The relation between quantum and thermal fluctuations has 
been studied previously via thermo field dynamics [8] under equilibrium and stationary 
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conditions. Here we treat the fully nonequilibrium problem. Using the microdynamics of 
a quantum system as starting point we view thermal fluctuations as statistical variations 
of the coarse-grained environmental variables with which the quantum system interacts, 
the exact microdynamics of the system and the environment obeying the laws of quan- 
tum mechanics. The problem under study can thus be stated equivalently as finding the 
uncertainty relation for an open quantum system [9]. 

On the third issue, in loose terms, one often identifies the high temperature regime 
of a system as the classical domain. On the one hand, one regards the regime when 
thermal fluctuations begin to surpass quantum fluctuations as the transition point from 
quantum to classical. On the other hand, from the wave picture (or consistent history 
formulation) of quantum mechanics we know that a necessary condition for a system to 
behave classically is that the interference terms in its wave function (or interference in the 
histories) have to diminish sufficiently and remain insignificant so that probability can be 
assigned to classical events [4,5] or that classical decoherent histories can be well-defined 
[6]. This is known as the decoherence process. Is there any relation between these two 
criteria of classicality? We show in this problem under the conditions studied that they are 
indeed equivalent: The time the quantum system decoheres is also the time when thermal 
fluctuation overtakes quantum fluctuations [see Eq. (22, 25)]. 

However we issue a warning here that this regime should not be called classical, as is 
customary in many quantum to classical transition studies. In fact, after the decoherence 
time only the first postulate of quantum statistical mechanics (QSM) is satisfied, the system 
can be described by non- equilibrium QSM. Only after the relaxation time, when the second 
postulate (more precisely, the implied version for canonical ensemble) is satisfied for the 
closed system, can one use equilibrium QSM. Classical has still a long way to go. It is well- 
known that quantum statistical effects can be important at very high temperatures (e.g., 
Fermi temperature for metals). This is due to exchange interactions of identical particles, 
a distinctly quantum effect. Only when the statistical properties of fermions and bosons 
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can be approximated by distinguishable particles, usually at high temperatures when the 
Fermi-Dirac or Bose-Einstein statistics approaches the Maxwell-Boltzmann statistics, can 
the system be rightfully called classical. In this regard quantum carries two meanings, one 
refers to the interference effect and the other refers to spin-statistics effect. 

To examine these issues, we use the simple Brownian model [4-6, 10-13] of a quantum 
open system and the influence functional method [11] to incorporate the statistical effect 
of the environment on the system. As the microdynamics is explicitly treated in this 
approach, one can study how the result depends on the properties of the bath and the 
system-bath interaction in detail. We consider an initial Gaussian wave packet (12) with 
nonzero averaged initial position and momentum and calculate the time dependence of 
the spread < Ax > 2 , < Ap > 2 and the uncertainty function Ut in the presence of both 
quantum and thermal fluctuations. The main result is given formally by Eq. (17), which, 
for an ohmic environment under the conditions of weak damping and high temperature, 
simplifies to Eqs. (19) and Eq. (24) respectively. As the initial state is arbitrarily- 
squeezed, we can show clearly how the squeeze parameter enters in these uncertainty 
relations [14]. 

Our system is a Brownian particle with mass M and natural frequency Qq. The 
environment is modeled by a set of n harmonic oscillators with mass m n and natural 
frequency uj n . The particle is coupled linearly to the nth oscillator with strength C n . The 
action of the combined system and environment is 

t r -, 

S[x,q}= fds ^-Mx 2 -^M02 x 2| + ^|l m ^2_ l m ^2 ? 2 I + J2{-C n xq n } (4) 

q L " n n 

where x and q n are the coordinates of the particle and the oscillators respectively, and fio 
is the (bare) frequency of the particle. We are interested in how the environment affects 
the system in some averaged way. The quantity containing this information is the reduced 
density matrix of the system p r (x, x') obtained from the full density operator of the system 
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and environment p(x, q; x'q') by tracing out the environmental degrees of freedom (q, q') 

+00 +00 

p r (x,x',t)= J dq J dq'p(x,q;x',q',t)S(q- q') (5) 



—00 —00 



The reduced density matrix evolves under the action of the propagator J r (x, x' ', t \ Xi, x[, 0) 
in the following way: 

+00 +00 

p r {x,x\t)= J dXi J dx\ J r (x,x',t I Xi,Xi,0) PrfaijX^O ) (6) 



-00 —00 



In general this is a very complicated expression since the evolution operator J r depends on 
the initial state. If we assume that at a given time t = the system and the environment 
are uncorrelated, i.e. that p(t = 0) = p s x p e , and that the environment is in thermal 
equilibrium at a temperature T, the evolution operator for the reduced density matrix has 
an analytic form given by [5,11-13]: (We use the notation of [5]) 



Jr(x f ,x' f ,t I x i ,x;,0) = Z (t)exp^{[u 1 (0)X i + u 2 (0)X / ]y i - [in(t) x i + Mt) x f\ Y f} 

-1 



x exp^jan^lf + [a 12 (t) +a 21 {t)]Y l Y f + a 22 (*)*/} 



(7) 

where X = (x + x')/2 and Y = x' — x. The elementary functions u a {s) are obtained as 
solutions of the following differential equations 



d 2 Uj(s) 
ds 2 

with boundary conditions 



s 

+ 2 J ds'r](s - s')ui(s') + fllui(s) = (8) 



ui(0) = l, ui(t) = 

u 2 (0) = 0, u 2 (t) = l 
and the ciij (t) are obtained from the integral 

t Sl 



(9) 



a>ij(t) = J ds! J ds 2 Vi(s 1 )u(s 1 - s 2 )vj(s 2 ) (10) 




where i?i(s) = v>2(t — s) and V2(s) = u\{t — s) . Here rj and z/ are the dissipation and noise 
kernels in the influence functional [11]. They are given respectively by 



+00 

/flUJ 
duj I(u) coth(— — ) cosujs (11a) 
2kl 



and 



+00 

v( s ) = ~r l( s )'i l( s ) = [ dui l^—L cosus. (Ho) 
as J uj 



c 2 

Here I(u) is the spectral density function of the environment, I(u) = ~ cu n ) 2m . 

n ™ n 

An environment is called ohmic I(u) ~ u, supra-ohmic I(u>) ~ u n ,n > 1 or sub-ohmic 
n < 1. The most studied ohmic case corresponds to an environment which induces a 
dissipative force linear in the velocity of the system. 

We now consider a Brownian oscillator with an initial wave function 

(x - x ) 2 



iK*, 0) = exp{- (X 4 J oJ + l - PQ x) (12) 

where cr is the initial spread, and xo,Po are the averaged initial position and momentum 
of the Gaussian wave packet. This gives p r (x,x',0) = ip^x, 0)ip(x', 0). One can calculate 
p r (xf,x'f,t) by performing the Gaussian integrals over Xi and x\ and get 

Pr(x f ,x' f ,t) =Z (t)exp{-^X T QX + B T x} (13) 

where the prefactor Zo(t) depends only on time. Here X = (Xf,Yf) T . The parameters 
xq,pq enter in the elements Bi of B but do not appear in the uncertainty function. The 
elements Qij(t),i,j = 1,2 of Q are given below. 

To calculate the averages of observables, it is convenient to use the Wigner function 
defined as 

W r (X,p,t) = J dYet pY p r (X-^,X+^,t), (14) 

The quantum statistical average of a system observable, e.g., x n or p n , with respect to the 
reduced density matrix p r {t) is given by 

< x n > T = J dxx n p r (x,x,t) = J dx J ^x n W r (X,p,t) (15) 
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where the subscript T indicates that the environment is at temperature T. Obviously the 
average contains both quantum and thermal contributions. We get the variances 



(Ax) 2 T =<x 2 > T -<x > 2 T = — 



; (Ap) 2 T =<p 2 > T -<p >r 



detQ 

Quit) 



(16) 



and arrive at the finite temperature uncertainty function 

h 



u T (t) = (A x f T (A P f T = n 2det ^ - *S9*M _ <M')<M'h 



[Qn(tW 

The matrix elements in (17) are given by 



iQn(t) Qu(t) Quit))' 



(17) 



Q22W 1 



+ 



2a< 



Quit) 4[u 2 (0)} 2 h[u 2 (0)} 2 

h 2 



[ui(t)] 2 au(t) - «i(0)«i^[a 12 (f) + a 21 (t)} 



+ 



4a 4 



+ W] 



Mo)] 



^— 2 {4au(t)a 22 (t) - [ai2(t) + a 2 i(t)] 2 } 



(18a) 



and 



Ql2(t) = ggig) 

Q11W QnW 



a 



iu 2 {t) f a 



&[u 2 (0)] 5 



^2(0) 

M 2 (*) 



Ml 



(0)«i(t)} 



(186) 



m 2 (0) I u 2 (0) 

The complete result can be obtained from solving the equations for the Ui{s) functions and 
the a,ij(t) coefficients. It is also instructional to look at some special cases where one can 
find simpler analytical expressions. 

For an ohmic environment 7(£) = 2-fo5(t) and ay(£), Ui(t) become harmonic and 
exponential functions. Let us introduce the following physical parameters: the effective 
frequency O = (fi 2 , — 7Q/4) 1 / 2 , the energy parameter e = = 1/t, the damping parame- 
ter a = ^ and the squeeze parameter 5 = 2n ^ a which measures the spread in the initial 
Gaussian wavepacket. For ohmic weak damping, a « 1, O ~ fio the uncertainty function 
is given for all temperatures to first order in a by 



^U T (t) ~ ^| e " 70 * +cothe 

+ coth e 
1-5 2 



1 - e" 70 * 



f (1 - S) 2 


1 - e" 70 * 


\ 45 





+ 



46 



sin 2Qt + a 



coth e — 



l-o- 
45 

1 + S 2 
26 



a sin20t \e~ lot 



(19) 



sin 2 Otl e" 27ot 



We see that there are two factors at play here: time and temperature. Time is measured in 
units of the relaxation time proportional to t re i = 7 C J~ 1 , and temperature is measured with 
reference to the ground state energy %VLq/2 of the system. There are also two parameters 
involved, a and 5. One can deduce various limits from this expression. For example, 
it is obvious that at t = 0, when the initial uncorrelated conditions is assumed valid, 
Ut(0) = h 2 /4, which is the Heisenberg relation (1). At very long time (t » 7 ( ^ 1 ), ^r(£) 
is insensitive to 5 and approaches U^ E as in (3) at finite temperature (as O ~ Q ). That 
means the Brownian particle approaches an equilibrium quantum statistical system. (For 
supraohmic bath this may not always be true). We can also see that for a T = bath 
(cothe = 1), £/r(t) has a leading term given by h 2 /4 (the Heisenberg relation) followed, 
for squeezed states S ^ 1, by terms of order a and a depicting both decay and oscillatory 
behavior. We understand that this is due to the action of quantum fluctuations alone. For 
a minimum- uncertainty initial state (5 = 1), we get for all finite temperatures 

h 2 f r 1 1 2 

U T (t) = —< e" 70 * + cothe |l - e" 70 *] \ +0(a 2 ) (20) 

Notice that there is no linear order damping term. 
At short times (t « 7 c j~ 1 ), 

U T (t) ~ U - [l + 2(5 coth A_ _ 1)lot + (t 2 )] (21) 

This simple expression is revealing in several aspects: The first term is the ubiquitous 
quantum fluctuation, the second term is the thermal contribution, which depends on the 
initial spread and increases with increasing dissipation and temperature. The time when 
thermal fluctuations overtake quantum fluctuations is (assuming the temperature is higher 
than the ground state energy) : 

ti = 1 —^a (22) 

2 7 o(5cothU-l) 

We expect this to be equal to the decoherence time scale tdec calculated for weak coupling 
for all temperatures in the study of quantum to classical transitions [15]. 
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At zero coupling, i.e., for an isolated harmonic oscillator, 

?2\2 



(23) 



This is the quantum uncertainty relation for squeezed states. The time-dependent term is 
the result of quantum dispersion. 

At high temperatures, To = 2 nci both the noise kernel and the dissipation 

kernel for ohmic dissipation become local, we have for the full range of damping strength, 

2 



h 



1 l rit) - r () 



1 -e 



-7o t 



+ T 



Q 2 
O 2 



(1-5) 
45 

1-5 2 
45 



2 r 



1 - e" 70 * 



sin 20t + a 



rl-<S 



a 



2 n 



T 



45 

1 + 5 2 
25 



sin20t - a 2 



T 



1 + 5 2 
25 



sin 2 Qt le _7ot 



sin 2 Ot I e" 270 * 



At short times (t << 7 1 ), 



1 + 2(r 5 - l) 7o t + 0(t 2 ) 



(24) 



(25) 



The time when thermal fluctuations overtake quantum fluctuations is given by 



(26) 



4kT-f 5 

This is equal to the decoherence time scale td ec calculated before for high temperature for 
the same model [4,5]. 

To complete the description of the evolution of the system, the second time scale of 
importance is the relaxation time scale, t re i = 7 - " 1 we referred to earlier, when the particle 
reaches equilibrium with the environment. It is at this time that equilibrium QSM can be 
applied to the description of the system. After this, for ohmic and subohmic environments 
the uncertainty relation takes on the Bose-Einstein form (3). At high temperatures the 
system reaches the Maxwell-Boltzmann limit and the uncertainty relation takes on the 
classical form (2). For supraohmic environments at low temperature, the highly nonlocal 

10 



frequency response may make it difficult for the system to settle down. The decoherence 
time scale is longer, and the relaxation can even be incomplete. This is the regime where 
one expects to find more intricate and interesting behavior in the interplay of quantum 
and thermal effects. The nonohmic results [16] and details of the present study [9] are to 
be presented elsewhere. 

After this paper was submitted for publication we learned of the work of Anderson and 
Halliwell [3] on an information-theoretical definition of uncertainty which, among reporting 
other interesting results, also confirmed our main findings here. 

We thank Juan Pablo Paz and Alpan Raval for discussions, and Jonathan Halliwell 
for explaining their information-theoretic results. This work is supported in part by the 
National Science Foundation under grant PHY 91-19726. 
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